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Abstract A dilute system of reacting particles transported by fluid flows is con- 
sidered. The particles react as A + A — > with a given rate when they are within 
a finite radius of interaction. The system is described in terms of the joint n-point 
number spatial density that it is shown to obey a hierarchy of transport equa- 
tions. An analytic solution is obtained in either the dilute or the long-time limit by 
using a Lagrangian approach where statistical averages are performed along non- 
reacting trajectories. In this limit, it is shown that the moments of the number of 
particles have an exponential decay rather than the algebraic prediction of standard 
mean-field approaches. The effective reaction rate is then related to Lagrangian 
pair statistics by a large-deviation principle. A phenomenological model is intro- 
duced to study the qualitative behavior of the effective rate as a function of the 
interaction length, the degree of chaoticity of the dynamics and the compressibil- 
ity of the carrier flow. Numerical simulations in a smooth, compressible, random 
delta-correlated-in-time Gaussian velocity field support the proposed heuristic ap- 
proach. 

Keywords Chemical reactions; dilute media; transport; large deviations; 
Kraichnan ensemble 



1 Introduction 

Many natural and industrial processes involve the reaction or collision of diffusing 
species transported by an outer flow. Such systems are typically modeled in terms 
of the reaction-diffusion-advection equation for the density field p. In the sim- 
ple case of the pair-annihilation reaction A + A — > 0, where two molecules react 
together to become inert, this kinetic equation takes the form 

d t p + V-(py) = -r P 2 + KV 2 p, (1.1) 

where K is the diffusion constant, v the velocity of the outer flow, and F the reac- 
tion rate. An important aspect withheld in such an approach is the assumed relation 
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between the microscopic stochastic rate, that is the probability that two given in- 
dividual molecules react, and the mesoscopic reaction propensity, that is the num- 
ber of reactions per unit time and volume written here as Fp 2 . To derive ( |1.1[ ) 
one must assume two properties to be satisfied at the coarse-graining scale from 
which the hydrodynamic limit is taken [30 j . First, volumes at the coarse- graining 
scale have to contain sufficiently many particles, in order to safely disregard finite- 
number fluctuations. Second, each particle within this volume must have an equal 
probability to react with all the others — well-mixing hypothesis. For this second 
condition to be satisfied, the coarse-graining scale has to be sufficiently small, i.e. 
smaller than the interaction distance between particle pairs. 

It is clear that in situations where the reactants are very dilute, the two condi- 
tions underlying might not be simultaneously satisfied. In particular, fluctu- 
ations due to a finite number of reactants might be so important to invalidate the 
mean-field assumptions leading to ( |1.1[ ). Much work has been devoted to model 
and study such fluctuations in situations where transport is negligible and diffu- 
sion dominates. It was shown that finite-number effects can be taken into account 
by adding an imaginary noise in the reaction-diffusion equation. The time evo- 
lution of statistical quantities is then obtained by averaging with respect to this 
noise (see II21II28I for reviews). Two equivalent approaches were independently 
developed using either the Poisson representation II 141 or a field-theoretical de- 
scription [9 24 J. In both cases, it is assumed that a reaction takes place with a 
given rate when two particles are on the same node of a lattice. Diffusion is then 
modeled by a random walk. After that, the hydrodynamic limit is obtained by 
performing the limit of vanishing lattice spacing. Such approaches, which were 
developed for diffusion-limited reactions, cannot be straightforwardly extended to 
cases where transport cannot be neglected. 

In many natural situations, the dynamics of very dilute reacting species is dom- 
inated by (possibly compressible) advection and their diffusion is almost negligi- 
ble. This is the case, e.g., of dust grains growing by accretion to form planets in the 
early solar system [19] or of phytoplankton confined in a two-dimensional ocean 
layers [25 1. Also, motile microorganisms in aquatic environments can detach from 
the fluid trajectories thanks to swimming and behave as tracers in a weakly- 
compressible flow [29]. Because of the large sizes of the transported species, 
molecular diffusion is negligible in all aforementioned examples. Another in- 
stance concerns the initiation of rain in warm clouds. Droplets, which are typi- 
cally occupying a volume fraction of the order of 10 -6 , grow by coalescence to 
form raindrops. Because of their finite size and mass, such droplets have inertia 
and their dynamics decouples from the flow. When they are sufficiently small, 
droplets behave as tracers in an effective compressible flow that depends on the 
local turbulence of the carrier airflow [22J. This leads to strong and violent fluc- 
tuations both in their local spatial concentration [27] and in the rate at which they 
collide lITTl . Classical approaches consist in using the (mean-field) Smoluchowski 
coagulation equation, which is a generalization of the kinetic equation to the 
case of a full population with a size distribution [26]. Many questions remain open 
on the statistical effects of turbulence on the timescales at which rain is formed [ 8 1 . 
In addition very little is known on the finite-number fluctuations induced by di- 
luteness that might also play an important role in assessing the average growth 
rate of droplets. The aim of this work is to develop a general framework which 
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can be used to address such issues. For that, we neglect diffusion and focus on 
particles transported by generic compressible flows. Also, we focus on the anni- 
hilation reaction A + A— >0 because we expect that this simple model will capture 
the essence of finite-number effects in advection-reaction systems. 



It is clear that, at long times, the particle number density decreases, so that 
this long-term asymptotics is very likely affected by finite-number effects. Kinetic 
approaches using < | 1 . 1 1 > predict that, ultimately, the average number of particles 
decays algebraically as f _1 Q. This behavior, which occurs at times much longer 
than the fluid velocity correlation time, can be explained in terms of an effective 
eddy diffusion. The flow compressibility might decrease the effective diffusiv- 
ity 1311 . but the latter remains anyway positive [17|. Hence, at long times, the 
spatial fluctuations of the concentration are smoothed out and a closed equation 
can be written for the spatial average of p. Here, we find that finite-number effects 
enhance the long-time decay of the average number of particles. For that we con- 
sider discrete particles that are tracers of the compressible fluid flow and that anni- 
hilate with a rate fi when they are separated by less than an interaction distance a. 
We show that the average number of particles does not decrease algebraically but 
rather exponentially as exp( — yt). This law pertains to the statistics of the relative 
motion between two reactants and the exponential decay rate y depends on both 
the microscopic rate fi and the flow statistical properties in a non-trivial manner. 
The main result of this work is the introduction of a novel Lagrangian approach 
in terms of non-reacting particle trajectories. We exploit these ideas to express y 
using a large-deviation principle for the time that two tracers spend at a distance 
below the radius of interaction a. 



The paper is organized as follows. Section [2]contains basic definitions and set 
the general framework of this work. The «-point number- density field is intro- 
duced and is shown to obey a hierarchy of transport equations. When integrated 
over space, these fields correspond to the factorial moments of the number of par- 
ticles that are present in the system. We then briefly discuss the zero-dimensional 
case and its relationship with standard studies of finite-number effects in well- 
mixed settings. In section[3] the hierarchy for the n-point number-density is solved 
by using a Lagrangian approach that consists in following the flow characteristics. 
At long times the particle number moments are shown to decay exponentially with 
a rate y that does not depend on their order. An analytic expression for y is given 
in terms of the Lagrangian statistics of (non-reacting) tracer trajectories. In section 
|4]we focus on the case of two particles. The exponential decay rate is then related 
to the time spent by the pair at a distance less than the interaction radius and is 
expressed in terms of the large fluctuations of the latter. Asymptotic arguments 
are then developed to relate the decay rate y to the microscopic rate fi . Section [5] 
contains an application to very dilute suspensions in smooth compressible Gaus- 
sian delta-correlated-in-time flows (belonging to the so-called Kraichnan ensem- 
ble; see, e.g., 1121 ). Finally, section |6]presents some concluding remarks and open 
questions. 
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2 Master equations in the presence of ad vection 

2.1 Settings 

We consider a set of particles {A ( } 1<KA ,(,) whose positions X, obey the equations 

^=v(X ; ,f), (2.1) 

where v is a prescribed li-dimensional differentiable velocity field with given 
isotropic, stationary, and homogeneous statistics. The velocity field v is assumed 
to have a finite correlation time and might be compressible with compressibility 



tr ( Vv T Vv) 

where the over line designates the average with respect to the velocity field re- 
alizations, and p 6 [0 : 1]. The extremal values p = and p = 1 correspond to 
incompressible and purely potential velocity fields, respectively. Whenever p > 0, 
the dynamical system defined by ( |2.1[ l is dissipative. Also, to maintain sufficiently 
generic settings, the dynamics is assumed ergodic and chaotic. Typically, when 
the flow is defined on a compact set, the trajectories generated by ( |2.1| l will con- 
centrate on a dynamically evolving strange attractor, while in the incompressible 
limit they remain uniformly distributed (see, e.g., lTT2l for more details on particle 
transport in compressible and incompressible flows). For the sake of simplicity, 
in this paper we will focus on the dynamics of tracers given by ( |2. 1 1 >. However, 
the results can be straightforwardly extended to a general dynamics given by the 
Newton equation X = F(X,X,f) as that ruling, for instance, the evolution of iner- 
tial particles ll22ll . 

The dynamics ( |2.1[ ) is supplemented by binary reactions between the particles. 
When |X, — X ; | < a the particles labeled i and j might react and annihilate (or be- 
come inert) . This happens with a rate jj. . As a consequence, the number of particles 
N(t) in the domain will decrease from its initial value N(0) = No in the course of 
time. We suppose that the radius of interaction a is smaller than the scale at which 
v varies. For example, in turbulent flows, this corresponds to assuming that a is 
smaller than the Kolmogorov length scale Moreover, in writing ( |2.1[ ) we have 
neglected particle diffusion. This is justified whenever the interaction distance a is 
larger than the Batchelor length scale £b — £k/VSc, where Sc = v/jc designates 
the Schmidt number — that is the ratio between the fluid kinematic viscosity V 
and the particle diffusivity K. For applications, our settings thus reduce to very 
large values of the Schmidt number. We notice that the Schmidt number can be of 
the order of thousands or higher in organic mixtures, biological fluids and gener- 
ically for particulate matter. Estimates based on Stokes-Einstein relation lead for 
a micro-meter sized spherical particle Sc w 1000 in air and 10 6 in water. 

Before discussing the master equations which describe the system, it is worth 
underlining that our settings involve different sources of stochasticity. First, the 
dynamics ( |2.1[ ) has to be supplemented by initial conditions on particle positions 
that we choose randomly. A second source of randomness comes from the intrin- 
sic stochasticity of the reaction process. Finally, the fluid flow is stochastic with 
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prescribed statistics. Consequently, ensemble averages corresponding to these dif- 
ferent sources of randomness must be considered; the corresponding notations 
will be introduced when needed. 



2.2 Master equations for ra-point number densities 

To fully describe all interactions between the particles in the considered system, 
all relative positions are simultaneously needed. It is then natural to define for 
1 < n < Nq (Nq = N(0) being the initial number of particles) the joint ra-point 
number density 

& n {x 1 ,...,x m t) = l £ 8(X h (t)-x 1 )---8(X in (t)-x„)\ , (2.3) 

where the brackets (•}» stands for the ensemble average with respect to reactions 
only and, by convention, the sum in ( |2.3| l is zero when N(t) < n. The particle 
indices vary between 1 and the number N(t) of particles in each realization of 
the reactions. The sum is not ordered so that, by construction, the spatial averages 
of the #«'s are the factorial moments of the total number of particles, namely 

Fi(f) = (n,/) d'jci = (#(/)}„, (2.4) 

F 2 (0 = JI^2(xi,X2,t)d d x l d d x 2 = (N(t)(N(t)-l)) fl (2.5) 



N(t) 



F»(f) = /•••/ ^n(xi,---,x n ,t)d d x l ---d d x n = ^ N ^^ n y J . (2.6) 

Hence, & n (t) jn\ is the average number of n-uplets among the N(t) particles still 
present at time t. 

Taking into account the transport by the flow v and the reactions, the joint 
ra-point densities (|2.3|l satisfy the following hierarchy of equations 



(9 & " 



ju£ / e(a-\xi-x n+l \)^ n+1 d d x n +i, (2.7) 
i=\ J 



where V X/ • is the ^-dimensional divergence with respect to the spatial variable x, 
and 6 the Heaviside function. The transport term comes from the identity 

^S(X ik (t)-x k ) = -V Xi - [8(X ik (t)-x k )v(x k ,t)] . (2.8) 

The two terms in the right-hand side are sinks accounting for the reactions: the 
first includes all binary reactions between any close-enough pair among the n 
considered particles; the second counts the reactions of any of these n particles 
with another one at a distance smaller than a. 
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The hierarchy ( |2.7| l for the «-point joint number density is exact and contains 
all information about the dynamics and statistics of the system. Moreover, it is 
closed at the order n = No and can be integrated by the method of characteristics 
(see Sect.|3]l. In the next subsection we relate this description in terms of factorial 
moments to standard descriptions used in the statistical physics of well-mixed 
chemical reactions. 



2.3 The zero-dimensional case 

When the hypothesis of well-mixedness is satisfied, as when for instance the inter- 
action radius is larger than the system size, transport does not play any role in the 
temporal evolution of the factorial moments \23\ . In this zero-dimensional 
limit, the hierarchy of transport equations \2.1\ becomes a simpler hierarchy of 
ordinary differential equations, i.e. 

= -^«(«-l)F„-|UnF„ + i, (2.9) 

which can be also derived from the standard master equation describing the re- 
action A + A — >0. Denoting with ^w(f) the probability of having N particles at 
time t, we have l30ll 

^ = | (N + 2) (N+ 1) 0> s+2 -^N(N- 1) *3V (2.10) 

By definition, the factorial moments are expressed in terms of 8?w{t) as 

M0=I^^v(0. (2.ii) 

Taking the time derivative of F„(f) and using ( |2.10| l yields to ( |2.9| l. 

Equation ( |2.1Q[ l is generally the starting point for describing the statistics 
of chemical reactions. In the limit of large number of particles, the mean-field 
appr oxima tion (N(t)(N(t) - 1)) M as (N(t) 2 )^ as (N(t))t is often used to write 
from( |2.10| a closed equation for the mean number of particles. This mean-field 
equation, obtained by replacing F2 by Fj in ( |2.9| l for n = 1, predicts a f _1 alge- 
braic decay of the number of particles. Clearly, the mean-field approximation will 
breakdown at long times when the number of particles is smaller. Many works 
have been devoted to describe fluctuations and deviations from mean field due 
to a finite number of particles [4 7 28 |. In particular, it is possible to show that 
solving equation \2.9\ is equivalent to computing the moments of the solution of 
a stochastic differential equation with a pure imaginary noise P1 [T5ll . However, 
this approach cannot be directly applied when the interaction radius is finite and 
transport is present. Reactions and spatial location of particles cannot be treated 
separately as it is usually done in reaction-diffusion systems, where the well-mixed 
hypothesis is assumed at the scale of the lattice spacing. The hierarchy of trans- 
port equations ( |2.7| l naturally couples the spatial distribution of particles and the 
reactions. The aim of the next section is to provide a formal solution of ( |2.7[ l in 
terms of Lagrangian trajectories (or tracers) of the carrier flow. These trajectories 
are completely defined by carrier flow and do not depend on the reactions. 
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3 A Lagrangian approach to particle reaction kinetics 

3.1 Finite-number closure and recursive solutions 

We now consider a finite interaction distance so that the presence of transport 
by the flow must be explicitly accounted for. Without loss of generality, we can 
assume that at t = the total number of particles is fixed and equal to No- As 
& n — for all n > No, the hierarchy ( |2.7| l is closed at the order No and can be 
integrated. 

Let us start with the joint A^o-point density 

^ + LVx,-[^v o v(x, : ,0] = -M^v o £ eia-lxt-Xjl), (3.1) 
m i=l i<j<N 

which can be formally solved in term of Lagrangian trajectories — that is using 
the method of characteristics. Denoting by Y(f;y) the solution to the differential 
equation 

^Y(*;y)=v(Y(*;y),0, Y(0;y)=y, (3.2) 
it is straightforward to check that the solution of ( |3.1| l is given by 
& No (xi,...,x No ,t) = 

r 1 N ° 

= /^V^ L ' <7 - Wo/ ° e(fl ~ |Y( ^^ 

■' i=l 

= (,^° e -^i<j<NoIo e ( a -\ Y ( s 'yi)- Y ( s -yi)\) ds \ (3 3) 

wo' 



with ^ o (yi,...,yjv ) = ^v (yi,...,yjv ,0). We have introduced here the so- 
called Lagrangian average (-)n , which is an average over all ,/Vo-uplets of non- 
reacting trajectories such that their location at time t is Xi,. . . , Xn . Such an aver- 
age is often used in compressible transport [ 16 1. We can now insert the expression 
for j£/v hi the equation for J^v _i obtaining 



N -i JVo-l 



" "™ 1 '■ E'vft-i'M]=-rt-i E e(«-|n-»/l) 



;'=l i<j 

I r N a -\ 



where 



M (^ e- A "o(" £ 0( fl -|Y( f ;y,-)-Y(?;y Wo )|)) , (3.4) 

\ ' =1 /No 

A n (t)=H £ f 0(a-\Y(s; yi )-Y( S ;yj)\)ds. (3.5) 

The solution of Eq.( |3.4[ l reads 

^v _ 1 (x 1 ,...,x iVo ,f) = (^ -ie- Aw o->«\ i _ (3.6) 



Wo-i 



'o-'W £ / 0(«-|Y(.;y ( -)-Y(5; yi v o )|)e^o->W-^oWd, 

< =1 Jo U 
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The general solution for ^ n can be obtained reintroducing ( |3.6| l in \2.1) solv- 
ing for and iterating until the order n. 

In writing ( |3.2| l we have implicitly assumed the existence and uniqueness of 
the Lagrangian trajectories and each trajectory is labeled by its initial position. 
This assumption is guaranteed as we have considered a smooth velocity field. For 
a fully turbulent flow and if a is much larger than the Kolmogorov length scale, the 
velocity field is no longer smooth but only Holder continuous with an exponent 
smaller than one. In such a case the solutions to ( |3.2[ ) are not unique and the 
deterministic Lagrangian flow breaks down. However, a statistical description of 
Lagrangian trajectories is still possible 1161 . 

We conclude this subsection by noticing that this new approach in terms of La- 
grangian trajectories allows us to directly apply known results from passive (non- 
reacting) transport (see [12] for a review) to address the problem of advection- 
reaction. These theoretical considerations will be mostly used in Sect. [4] 



3.2 Long times and very dilute closures 

When considering very dilute systems and small interaction radii, the hierarchy of 
transport equations can be closed at each order as the second term of the right hand 
side of Eq. \2.1\ can be neglected with respect to the first one. In particular, the 
closure becomes exact in the long time limit when the mean number of particles 
depends only on the statistics of the relative distance of one pair of Lagrangian 
trajectories. 

The long time behavior of the joint «-point density can be directly obtained 
from equation \3.6) . Noticing that A^ _ i (s) — Apj (s) < 0, at the leading order, the 
joint n-point density is given for n > 2 by 

& n (x h ...,x n ,t) = (c- A »V) n (3.7) 

Integrating over all spatial variables, for n>2, we obtain 

/ „ -M EJge(a-|Yi«-Y;(i)|)dA 

F„(0 = (^„°e "J j (3.8) 

where now the average is over all possible n-uplets of tracers. Finally, averaging 
the factorial moments over the realizations of the velocity field, we can define the 
effective reaction rate as 

% = -lim-lnML (3.9) 
t^t F„(0) 

Note that the average over different realizations of the carrier flow, denoted with 
an overline, is performed before taking the logarithm. This important point will be 
further discussed in Sect. 14. II 

In the zero-dimensional limit, by definition the Heaviside function inside equa- 
tion ( |3.8| l is equal to one, and thus the effective rate is simply given by y n = 
\in(n— l)/2, as confirmed by numerical data (not shown). When the system is 
not well-mixed at all scales, the effective rate (|3.9|l cannot be computed for generic 
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Fig. 3.1 Temporal evolution of (N p ) for p = 1,2,3,4 and initial values Nq = 2 (black solid 
lines) and No = 100 (red dashed lines). Simulations are obtained by Monte-Carlo simulations 
of the reactions using a delta-correlated in time velocity field (see Sec. [5] for settings of the 
flow). Average is taken over reactions, realizations of the flow and initial positions of particles. 
The parameters of the simulation are fi = 1, p = 0.1 and a/L = 0.015. The red dashed lines 
have been arbitrarily multiplied by a factor 1000 for the sake of representation. Inset: temporal 
evolution of the inverse of the average number of particles for ^0=100. It displays the mean-field 
prediction (Nty))^ ~ JVb/( 1 + CNq t) which is indeed expected to work at short times. 



flows as multi-time and space correlations of Lagrangian trajectories are involved. 
However, an important role is played by the case n = 2 that only depends on the 
two-point motion. Indeed, it is easy to check from equation ( |3.8| l that y n is an in- 
creasing function of n. In addition, the moments of the number of particles can 
be expressed as a linear combination of the factorial moments. The rate 72 is thus 
always the leading exponent and for all p > 1 we have 

(N p )~e-K<. (3.10) 

This exponential decrease is expected to be valid at long times regardless of the 



initial number No of particles considered. This behavior is evidenced in Fig. 3.1 
where the temporal evolution of the moments (N p ) is displayed for No = 2, 100 
and different orders p. For No = 100, the long-time departure from the <=c t~ x mean- 
field prediction is clear from the inset. Independently of the order of the moment 
considered, the long-time dynamics depends only on the statistics of the relative 
distance of one pair of tracers, which motivates the next section that is devoted to 
the study of the two-point motion. 



4 Two-point motion 

4.1 Effective rate 

As the long-time behavior of the system is dominated by the two-particle dynam- 
ics, it is natural to close the hierarchy \2.1\ at the second order and to look at the 
mean numbers of particles and of pairs: 



n 1 (t) = (N(t)) ll = F 1 (t), »2(0^M0M')-1]> M = ^F 2 (/). (4.1) 
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The overbars designate the average with respect to the realizations of the velocity 
field v, which is assumed statistically homogeneous in space so that the overbars 
include the average over the particle initial positions. Both quantities, n\(t) and 
n%(i), are expected to exponentially decrease with the same rate "fi, which here 
and in the following will be denoted by y. An effective equation describing such 
a behavior is simply given by 

d«i d« 2 , . r,. 

~d7 = ~ 7 " 2 ' ~d7 = ~ yH2 ' ^ 

The effective rate y depends on the microscopic rate /i, the interaction radius a 
and the statistics of the carrier flow trough equations ( |3.8[ ) and ( |3.9[ ). However, 
these equations do not allow us to readly obtain an explicit expression for y. The 
dependence of such an expression on the parameters of the system can be obtained 
by using standard tools of dynamical system and statistical physics. 

The effective rate y, defined by equation fl3.9) , depends on the particles tra- 
jectories only through and thus is determined by the two-particle statistics. 
Consequently, without loss of generality, we will focus on a system initially hav- 
ing Nq = 2 particles. For such a system, = 2 and the average number of pairs 

reads 

n 2 {t) = e -M/o , e(a-KW)d5 j (43) 

where R(t) = |R(?)| = |Yi(f) — Y2(?)| and the average is over all realizations of 
the carrier velocity field v. The argument of the exponential in Eq. ( |4.3| l contains 
a time integral that, for ergodic two-point dynamics, self-averages, i.e. 

1 r* 

lim-/ G(a-R{s))ds = P^{a)=F{R(s) <a\, (4.4) 
f Jo 

where P is a shorthand notation for probability. This observation leads to the fol- 
lowing naive prediction 

ynaive = M^2 < ( fl )- ( 4 - 5 ) 

This estimate is consistent with what could be expected on a heuristic basis, how- 
ever it is valid only when considering a fixed realization of the carrier flow. This 
kind of setting corresponds to the so-called quenched disorder in statistical me- 
chanics. Indeed, let us fix a realization v of the carrier flow and the initial positions 
Xi(0) and X2(0) of the pair of reacting particles. If we now denote by N(t) the 
variable which is equal to 2 when the two particles have not reacted and other- 
wise, we clearly have that 



lim I rn(iV(0 [N(t) -!]>„=- Hm £ jf 6(a - R(s))ds = -/i B 



Notice that, here, the average is performed solely with respect to reactions. The 
quenched effective rate is thus given by \A.5\ . The average over fluid realiza- 
tions can be subsequently performed but it does not modify the result (j43J as 
]n(N(t)[N(t) — l])n is a self- averaging quantity. 

The situation is quite different when considering an annealed disorder, that 
is when the fluctuations of the carrier flow are taken into account to compute the 
number of pairs, according to the definition ( |4.3[ ). In this case, as ni is not a self- 
averaging quantity, fluctuations of the carrier flow affect the long-time behavior 
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and, hence, y is not simply given by ( |4.5[ ). In most physical situations, fluctuations 
of the carrier flow and randomness of the initial particle positions cannot be dis- 
regarded, so that annealed averages must be performed to define the effective rate 
( |3.9[ ). As discussed in the next subsections, in this case, the effective rate yean be 
related to the large deviations of the fraction of time spent by two particles within 
a distance less than a. 



4.2 Large deviations of the reaction time 

Determining the long-time asymptotics of «2 requires to understand the behavior 
of the finite-time average 

© = - f 6(a-R(s))ds. (4.6) 
t Jo 

Under ergodicity hypotheses on the two-point dynamics, this self-averaging quan- 
tity converges at long times to P^~(a). When the dynamics is sufficiently mixing, 
the deviations of from its average are described by a large-deviation princi- 
ple 1101 . The probability density function (PDF) p(&) thus has the asymptotic 
behavior 

-]imj\np{©) = Jf(0), (4.7) 

where is a convex positive rate function, which attains its minimum, equal to 
zero, at — Pf (a). Plugging this form in ( |4.3| ) we obtain that 

n 2 (f)« r e -^ 0+ ^ »d0. (4.8) 
Jo 

A saddle -point estimate of the integral shows that the effective rate, 

y= -lim-lnn 2 (0 = mf \fi0+Jtr{0)], (4.9) 

is the Legendre transform of the rate function M The convexity of M" implies 
that y is a non-decreasing function of fi . 

When the reaction rate fi is small compared to the scale of variation of Jff , the 
infimum in \A.9\ is attained very close to the minimum of M 3 at = (a). There, 
as a result of the central limit theorem (CLT), the rate function is approximately 
quadratic J4?(®) ~ (0 -P 2 < (a)) 2 /(2a 2 ) and the infimum in ( f4~9] ) is attained for 
0=0*= P 2 < (a) - jua 2 , so that 

y~;UP 2 <(a)-;U 2 y, (4-10) 

where a 2 /t is the variance of the fraction of time spent by the pair at a distance 
less than a. The expression ( |4.10[ ) is valid as long as 0* is close enough to P^i®) 
to neglect the sub-leading terms in the quadratic approximation for M 1 , i.e. in the 
domain of validity of the CLT. Notice that ( |4.10| ) requires \1<^P^ (a) /a 2 . Observe 
also that the naive estimate ( |4.5| l for a quenched velocity field only captures the 
linear behavior of y close to fi = 0. 
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The asymptotic behavior of y at large values of fi relates to that of the rate 
function ffl at small values of the fraction of time spent below a, indeed, from 
( |4.9| ) we have 

lim y = lim JT(0). (4.11) 

The events leading to small 0's correspond to those in which the particle-pair 
separation remains larger than a for a very long time. More precisely, the value 
of ffl at = relates to the distribution of the random time Thit needed by pairs 
that are initially far apart (e.g. at separations of the order of the system size) to 
approach each other at a distance less than a. We can indeed write 

P(0 < e) = P(0 < £ | T m > t)W(T m > t) +P(0 < e \ T hit < t) P(T Wt < r). (4.12) 

When r^t > t, the two particles have never been at a distance less than a, which 
implies P(0 < £ | T^it > = 1 f° r an Y £ > 0. Conversely, when T^n < t, the fraction 
of time is finite and P(0 < £ | r hit < t) -> for £ -> 0. Hence, it follows that 

J^(0) = -limylnP(7h it >f). (4.13) 

Different behaviors of the effective rate are thus expected depending on the tail 
of the distribution of 7hi t . For the proposed large-deviation approach to be valid, 
one expects P(7hit > t) to decrease at least as fast as an exponential. Otherwise, 
one would obtain Jrif (0) = 0, violating the convexity of 3ft 1 . When P(7h; t > f) de- 
creases faster than an exponential, J^f(0) = °°, and the effective rate y — > °° for 
ji — > oo. The way it diverges cannot be obtained from the distribution of only, as 
it depends on the functional form of Jj? close to = 0. In the intermediate case, 
when P(7hit > t) decreases as an exponential, 3f (0) is finite and the effective 
rate approaches a finite value "fa = Jff(0) when fj. — > oo. Depending on whether 
H* = — dJrf? /d0\@=o is finite or not, / saturates to for /I > /I* or approaches it 
asymptotically. Note that an exponential behavior for the tail of the distribution of 
Thit is not a particular case but is actually expected to be a quite generic circum- 
stance. Indeed, when considering times much longer than all relevant timescales 
of the flow, memory loss implies that the hittings of a form a Poisson process, and 
the hitting time is an exponential random variable. 

In the next subsection we introduce and study a phenomenological model for 
which the rate function can be explicitly computed. 



4.3 A phenomenological model for reactions in bounded chaotic flows 

The essence of chaotic flows is the competition between stretching and folding. 
Stretching results from the fact that close trajectories separate exponentially at a 
rate given by the largest Lyapunov exponent A . Folding, which is generally due to 
boundary conditions, prevents the inter-particle distance to grow indefinitely and 
is necessary for the system to converge to a statistical steady state. Qualitatively, 
the time evolution of the distance R between two particles resembles the trajectory 



shown in Fig. 4. 1 The particle pair spends long times at a distance of the order of 



the size L of the domain and performs rare excursions to very small separations. 
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Fig. 4.1 Illustration of the decomposition in hitting events of duration 7j|j t and escapes of dura- 



When interested in particles reacting within a distance a -C L, it is quite natural 
to decompose the inter-particle dynamics in a sequence of consecutive phases, 
each corresponding to a given approaching event and separated by returns to L. 
Starting from R = L, a first time subinterval is required to hit a for the first time. 
Then the particle separation grows to touch again L. The time interval [0,/] is 
thus decomposed in n rea random subintervals of length TjJL + 7" e "c' where 72, is the 
hitting time from R = L to R = a and r e " c the time needed to reach L for the first 



time starting from a (see figure 4.1 1. Because of the memory loss occurring when 



the pairs are at a distance L, these sub-intervals can be considered independent 
and identically distributed. The first substage corresponds to many recurrences to 
R = L during time lengths that are independent from each other because of the 
imposed boundary condition. The hitting time is dominated by the sum of these 
recurrence times and this leads to assume that TjJL has an exponential distribution 
of mean 7hit- Concerning the escape time T" sc , it is clear that when a <C L, its 
behavior is dominated by the long-time convergence to the Lyapunov exponent. 
We can thus assume that it is deterministic, namely r e " c w r esc = (1 /A)ln(L/a). 
In addition, in this limit, we always have r esc <C r h " t and thus 



T" 1 



hit - 



(4.14) 



According to the discussion of the previous subsection, the times 7]jJ t are exponen- 
tially distributed and « rea is a Poisson random variable of mean f/Tirit. 

During the n-th interval of length r e " c , the particles spend a time a„ T" sc below 
a distance a. The total fraction of time spent at R < a can thus be written as 



1 n rea 
1 n=l 



i —ST 



(4.15) 



where a designates the average fraction of time spent below a during r esc . We 
have assumed in this model that the fluctuations of a do not play an important 
role. Notice that under these assumptions, P^(a) = = aT esc /Th\ t - We see from 
equation ( |4.15| l that is proportional to n rea . Its probability density is thus given 
by the Poisson distribution. Using the Stirling formula when t — > °°, we obtain 
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It is easy to check that Jff(©) given by (|4.16[ ) is a convex function that vanishes at 
its minimum reached at @ m i n = oT e sc/7hit = P^(a), as expected. Also, one sees 
that Jf?(0) = l/7hit and Ji?'(Q) = — °°, where the prime denotes the derivative 
with respect to ©. In the light of the considerations of the previous section, this 
implies that the effective reaction rate y asymptotically converges to y„ = 1 /T^t 
when fi — > oo. With the explicit form ( |4.16| l for the rate function we actually obtain 



1 

Thit 



1 



(4.17) 



Note that this expression also reproduces the small-/i asymptotics ( |4.10| l with 

(T 2 =P<(a) 2 Wf 

The explicit form ( |4.17| l relies on the strong simplifications of the particle dy- 



namics that we have made (cf. Fig. 4.1 1. However, this model is particularly useful 
to understand how y depends on the reaction radius a, the flow compressibility, and 
its chaoticity that will affect both and 7hi t . The time 7hi t + r esc ~ 7h; t can be 
seen as the recurrence time from a to a. Therefore, by the Kac Recurrence Lemma 
0201 . it is expected that T^it ~ l/^2 < ( a )- However, the proportionality factor de- 
pends on the Lyapunov exponent and on the compressibility. More precisely, by 
definition we have that 

7^P 2 < (fl) = ar esc = |ln(L/fl), (4.18) 

so that in this model and Thu can be explicitly related if a is known. 

Determining a corresponds to computing the average fraction of time spent 
by a particle pair at separation R(t) < a during a time interval r esc given that the 
initial is separation R(0) = a. We can estimate this quantity by approximating 
p(t) — ]n(R(t)/a) as a Brownian motion with drift A and noise variance A — 
which typically well describes the long-time dynamics of the relative separation 
between particles advected by smooth finite-correlation-time flows 121 1121 — , and 
thus computing the fraction of time spent in the negative half line, p(t) < 0, during 
a time interval r esc . Exploiting known results for Brownian motion with drift — 
namely, the generalized arcsine law HI (see Appendix [Ajl — we can obtain an 
analytical expression for the PDF of a, from which we derive that 

a = <P ((X 2 /A)T esc ) = <P {{X/A)\n(L/a)) . (4.19) 

The function cannot be easily expressed analytically but can be numerically 
computed with arbitrary precision. It strictly decreases from 1/2 to zero at infin- 
ity. In addition, the assumption on the relative particle dynamics implies that the 
separation R at time f between two particles which start at time with a separation 
Ro is distributed according to the lognormal distribution 



mt) = l Rvkl^ 



[ln(R/R )-Xt] 
2At 



(4.20) 



In particular, equation ( 4.20 1 implies for the generalized Lyapunov exponent L(q) = 

hm,_ i . {=0 lim So _ > .o(lA) m {R/Ro) q = +q 2 A/2 |6)- For a smooth velocity field, if 
the dynamics is sufficiently mixing, statistically reversible and stationary, there 



15 




Fig. 4.2 A) Average hitting time Ty, normalized by P^) for L/a = 2000 and as a function 
of the solid line correspond to the prediction \4.2J>\ using the generalized arcsine law. B) 
Dependence of the effective rate y (normalized by X ) on the correlation dimension &i f° r A 1 / X = 
land L/a = 2000. The solid line is obtained by using equation \4.17\ and 7hjt given by equation 
(4.23) . In both figures the symbols show data obtained from numerical simulations using the 
Kraichna n flow (see sectionpl. The constant c = 1.125 comes from the proportionality factor in 
equation j4.23fr . 



is a general result stating that for a certain exponent £$2 > 0, (R(t)/Ro)~^ 2 is a 
conserved quantity and thus L(— £^2) vanishes [18,3]. The exponent £^2 coincides 
with the correlation dimension — a classical observable in dissipative dynami- 
cal systems (see, e.g., 11231 ) — that rules the scaling dependence of the stationary 
probability of finding two particles at a distance less than a 

P 2 < (a)~ {a/Vf 2 for a<L. (4.21) 

For the lognormal distribution ( |4.20| l we have that L(— 2k / A) = 0, so that 

9 2 = 2X/A. (4.22) 

Using such a relation, the average fraction of time a ( |4.19| l can be re-expressed 
in terms of S>2- These considerations, together with the relation ( |4.18| l, lead us to 
express the average hitting time as 

17^- (L/a)* In (L/a)* [(^/2) In {L/a)]. (4.23) 



As seen in figure 4.2 \, these approximations are consistent with numerical data 
obtained in random Gaussian compressible flows. Such flows are chosen within 
the Kraichnan velocity ensemble; we will come back on such numerical inves- 
tigations in the next section. Note that using the above expression for 7h; t and 
equation ( |4.17| l, the effective rate 7 normalized by the Lyapunov exponent X can 
be expressed as function of QI2, a/L and ju/A. The dependence of y/X on £^2 is 
shown in figure |4~2] B as a solid line and is again in good agreement with the nu- 
merical data in random flows, providing support to the approach described above. 
Note that there is a clear exponential decrease with 3%. In general, the correla- 
tion dimension depends on the flow compressibility p. For incompressible flows 
(p — 0), the trajectories are space-filling and ^2 = d, the space dimension. For 
compressible flows, Q>2 is a decreasing function of p. We thus have that, generi- 
cally, compressibility enhances reactions. 
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The model presented in this section, although phenomenological, is based on 
quite general assumptions, such as the existence of a finite correlation time and 
of at least one positive Lyapunov exponent. The precise shape of the rate function 
( |4.16[ ) can vary from one system to another when more realistic dynamics is con- 
sidered. It is however reasonable to expect that the general qualitative properties 
of the effective rate /remain essentially the same. 



5 Applications to very dilute suspensions in compressible Kraichnan flows 

To illustrate and further validate the approach developed in the previous section, 
we focus here on two-particle motion in a special class of random velocity fields, 
that is the smooth (compressible) Kraichnan ensemble [16 12]. The main advan- 
tage of such a class of stochastic flows is that, thanks to time uncorrelation, we 
have analytical control on most of the quantities we need. Further it is very effi- 
cient to implement numerically so to gain high quality statistics. 



5 . 1 The Kraichnan velocity ensemble 



We focus on dimensions d > 1 . The velocity v is a homogeneous, isotropic, Gaus- 
sian spatially smooth field, with zero mean and two-point correlation function 



v i (0,t)v j {r,t')=2D i j(r)S{t-t'), D, 7 (r) = A)%-rf l7 (r)/2. (5.1) 
r| — > the spatial part of the correlation is prescribed 



For small separations r 
according to 11161 



dij(r) =/>![(</ + !- 2p)5 u r 2 +2(pd - l) ri rj\ . 



(5.2) 



The flow parameters are as follows: the diffusivity Do controls the single parti- 
cle dispersion properties; Di has the dimension of an inverse time and sets the 
intensity of velocity gradients; finally, p £ [0 : 1] controls the compressibility de- 
gree as defined in Sect 2.1 its extreme values correspond to incompressible and 



potential velocity fields, respectively. When advecting particles with the velocity 
field defined by Eqs. ( |5.1[ )-( [5T2) l, we will always consider the presence of boundary 
conditions, which are necessary in order to ensure statistical steady state proper- 
ties (see also below the description of the numerical implementation). 

Particles evolving in such flows with the dynamics Eq. ( 2. 1 1 define a stochastic 



dynamical systems of which we know explicitly the Lyapunov exponents, A# = 
Di [d(d - 2k + 1) - 2p(d + {d — 2)k)] (k=l,...,d), and the rate function of the 
stretching rates, which is quadratic (see the review [12] for a summary of the 
known properties of Kraichnan flows). From the latter result it is possible to prove 
that the lognormal distribution ( 4.20[ > for two particle separation becomes exact 
at sufficiently long times with A~=A\ = D\(d — \ )(d — 4p) and A = 2D\(d — 
1)(1 +2p). In particular, for the Kraichnan model, Eq. (4.22 1 (see discussion in 
Sect. 14.3 b holds true with 



(5.3) 
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It is worth noticing that $>2 — d for p = 0, which simply means volume preserv- 
ing dynamics in the incompressible limit, and that ^2=0 for p c = d/4, which is 
the value of the compressibility p above which the first Lyapunov exponent be- 
comes negative. The latter property indicates that above a critical compressibility 
p c there is a transition to a regime of particles trapping with particles asymptot- 
ically collapsing onto a single point see, e.g., Ref. 11161 . Notice that this regime 
cannot be realized for d > 4. 

The intuitive content of the lognormal behavior ( |4.20 1 is that the log of the 



separation between the particles, In [/?(/")], essentially behaves as Brownian motion 
with drift A and diffusion coefficient A/ 2, but for the presence of boundaries at 
R — L. We will see in the following subsection that this is a very good approxima- 
tion, which will be exploited to model the two-point quantities we are interested 



in, similarly to what has been done in Sect. 4.3 

In the following we will focus on the two-dimensional (d = 2) compressible 
Kraichnan model, which is the simplest non-trivial example and is quite efficient 
for numerical implementations. In particular, we shall consider compressibility 
values p< p c = 1/2 to stay away from the trapping regime, which is trivial from 
the point of view of reaction dynamics. 

As for the numerical implementation we consider two possible schemes. The 
first, which is useful when considering many particles (this is that used for Kg. [3.1) , 
consists in generating a random velocity field, with periodic boundary conditions, 
as the sum of independent Fourier modes with |k| < 2 and amplitudes such that 
Eqs.( |5.1| )-( |5.2| l are satisfied at small scales. When interested in two-particle mo- 
tion, we use a second scheme, which is much more efficient. It evolves directly 
the stochastic equation for the separation R between the particles. We use a dis- 
cretization by means of a standard Euler-Ito scheme that reads 

Ri(t + dt) = Ri(t) + V2dtVi(R,t) for i=x,y 
with V x = y/Di(\+2p)R x r\ x - y/Di{3-2p)R y j] y 

V y = y/Di(\+2$>)R y T\ x + y/Di{3-2p)R x T} y , 

which only requires to generate two independent zero-mean Gaussian random 
variables, f]i. The above equations indeed prescribe that V,(R,f)V/(R,f) = dij(R), 
with djj given by Eq.(5.2l for d = 2. Finally, to ensure a statistically steady dy- 



namics, reflective boundary conditions are imposed at \Rj\ =L/2. In principle, this 
scheme can be generalized also to more than two particles 11131 . 



5.2 Numerical results 

As discussed at length in the previous sections, the long-time statistics of very 
dilute systems is dominated by the two-particle effective reaction rate y. This is 
confirmed by numerical simulations of the compressible Kraichnan flow: as shown 
in Fig. 3.1 the decay of the moments of the number of particles is indeed asymp- 
totically dominated by y. Therefore, here, we focus on the two-point dynamics 
in two-dimensional Kraichnan flows with the aim of illustrating numerically the 



heuristic picture set up in section 4.3 
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Thanks to the Lagrangian formulation of the problem (see Sect. [3] and also 
Eq. (4.3 l) we do not need to implement the reactions and it is enough to study the 
statistics of the fraction of time = (1/f) J 6(a — R(s))ds spent below the reac- 
tion distance a. We can start with considering the average of 0, which is nothing 
but the probability P^i®) to find the particle pair at separation below the reac- 

a/L)^ 



tion distance. For small a, we have that P^ (a) 



(see Eq. 



4.21 



with the 
the self- 



correlation dimension &2 given by (5.3 1. Concerning the fluctuations o 
averaging quantity 0, clearly its variance decays in time, i.e. ((0 — P^(a)) 2 ) = 
G 2 /t. Also, its probability distribution obeys at long times the large-deviation 



principle (4.7 1, i.e. p{0) 



one observes that a <x 



P< 



exp[—t,J4?(0)]. As shown in figure [5TT] for small a, 
(a), so that the scaling behaviors of both P^{a) and 
a 2 are controlled by the same quantity, namely the correlation dimension as 
well verified by the numerics. (Notice that consistency between ( |4.17| l and ( |4.10[ ) 
would require <7 2 = (a) 2 7hj t , which is P^ip) as ^Mt ~ 1 this ls ver ~ 

ified in Fig. 5.3 ) Moreover, figure 5.2 shows the rate function $f computed at 
different times for a flow with p = 0.1 (different values of the compressibility 
provide qualitatively similar results). The collapse of the curves demonstrates the 
validity of the large-deviation principle. It is also worth noticing that for tiny de- 
viations from the average the rate function is well approximated by the parabola 
(0 — /^(a)) 2 / (2a 2 ) associated to the central limit theorem. However, it is inter- 
esting to note that for larger deviations there is a clear departure from this behavior, 
and the rate function becomes asymmetric. This is consistent with the phenomeno- 
logical expression (4.16 1. 

In order to test the phenomenological picture discussed in section [43] we need 
to study the statistics of the hitting time T^u which is at the base of our reasoning. 
To this aim we perform a large number of numerical experiments where we set 
R(0) = L/2 and measure the first time T^ t such that Z?(7iut) = <*■■ The distribution 
of such times is shown in Fig. 5.3 for p = 0.1 and several values of a. We expect 




a/L 

Fig. 5.1 (Color online) P 2 < ( a ) (°P en symbols) and a 2 (full symbols) as a function of a/L for 
different compressibility value s as labeled. The solid lines display the power law (a/L)® 2 '^, 
with &2{P) given by Eq. {53}. 
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Fig. 5.2 (Color online) The rate function J>P as a function of © fora/L = 5- 10 3 and p= 0.1 at 
different times as labeled. The solid black line denotes the central limit theorem approximation 
Jif(0) = (0 — P^(a)) 2 /(2a 2 ). The arrows show where the minimum Jf? = is attained, i.e. 
for O = P^(a), and the region of validity of the approximation. 



such a distribution to be exponential, P(7h; t ) ~ exp(— 3iit/^hit)/^hit- As shown 
in the figure, this is rather well verified. Nevertheless, deviations are observed 
when a/L is not too small. Indeed, an exponential PDF for corresponds to 
assuming that the hitting events follow a Poisson distribution, requiring the hitting 
process to be memoryless. Co mplete memory loss can be safely assumed when 
Thit 3> 1/A which, in virtue of (4.18 1, means 7hi t A = ln(L / 'a) / P^ (a) <C 1 so that 
the larger the compressibility the smaller should be a/L for the exponential PDF 
to be a good approximation. As shown in the inset of figure |5.3| we indeed have 



that A7hit ~ 1 /P 2 < (o)> providing a further support to ( |4.18[ ) and (4.23 I. This is also 



consistent with the results shown in figure 4.2 



In principle, we can now plug the measured T^n and P^ (a) into equation (4.16 1 
and compare the phenomenological prediction with the measured rate function 
Jtf '. Although qualitatively capturing the main features of Jtf?, the approxima- 
tion does not really match the numerical data (not shown). Notwithstanding, the 
measured rate y is well captured by the formula (4.17 1 obtained from our heuris- 
tic arguments. Clearly, this approach catches the qualitative behavior of y and as 
seen in figure |5.4| it is also quantitatively correct both in the limit /I — > where 



5.4 3 



y~ ixPjia), and in the limit ji — > °° where y« 1 /Thit- In particular, figure 
and its inset show how well the prediction ( 4. 10| ) coming from the central-limit 
theorem works for small values of ji . 

In section 4.3 we have seen that approximating the evolution of ln[R(t)/a)] as 
a Brownian motion with drift leads to a predictio n for the normalized rate y/A 
in agreement with numerical data for fi = A (Fig. 4.2 1. This suggests the possi- 
bility of exploiting similar ideas in order to improve the approximation for y at 
intermediate values of jj.. 
As sketched in Fig. 



4.1 we decompose a time interval t in n lea intervals of 



length r h " t + r e " c . The hitting times from L to a, T^L are random, exponentially 
distributed variables (see Fig. 5.3 I. Moreover, we can assume that at time t the 
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2 4 6 8 _10 12 14 16 18 

T hi /Thi, 



Fig. 5.3 (Color online) PDF of the hitting time normalized with the average hitting time T^ti for 
p = 0.l for several reaction distances a, as labeled. Notice that for a/L < 5 • 10~ 3 all curves 
collaspe onto the exponential PDF (in blue), deviations from pure exponential PDF can be seen 
only for larger values of a/L (in red), which however possess a clean exponential tail. Inset: nor- 
malized average hitting time A7],it vs a /L (symbols) comparared to the inverse of the probability 
P 2 < (a) . Notice that for small-enough distances 7^ <* 1 /P£ (a) ~ (a /L)- @1 W . 




Fig. 5.4 (A) 7 vs fi for p = 0. 1 and a/L = 5 • 10~ 3 the solid lines shows the small and large fi 
asy mptotic s, respectively. (B) y/fl vs fi the solid line displays the central-limit theorem predic- 
tion l |4.10} , which is made more clear in the inset by eliminating the constant term. 



distribution of n rea is Poissonian, P(n rea ,t) = K""* exp(— A r )/n rea !, with mean K — 
f/7hi t . Let us denote by a" the fraction of r e " c during which the particle pair is 
within a distance a. We can write / ' 9 (a - R(s))ds = £j^ a a n T^ sc so that 



e-f= e -MJi|8(«-*W)= £ P(«aa,0 

n rea >0 

= £ P(n iea ,t)[e-^]" Ka = £ P(« rea ,0/3"-=e-( 1 -/ 3 ^=e"'^. 

«rea>0 «rea>0 

In the above derivation we made essentially no approximation but the statistical 
independence of the time subintervals. This assumption is justified by the memory 
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loss of the process. As a result we have an explicit expression for the effective 
reaction rate 

y=(i-j8)/5S=y..(l-P), (5-4) 

which is consistent with the asymptotic behaviors for small and large ju. Now we 
need to compute J3, which can be expressed as 

P = e -M«r esc = r dT esc P(T esc ) [ d(ar' lor «P(a|r MC ) , (5.5) 
Jo Jo 

where P(T esc ) is the PDF of r esc and P(a|r esc ) the PDF of the fraction a con- 
ditioned on r esc . Next, as in section |43| we assume r esc to be non-random with 
P(T esc ) = 5(r esc — A -1 ln(L/a)). Also, exploiting again the approximation of the 
dynamics of ln(R(s)/a) as a Brownian motion with drift, we use t he generalized 
arcsine law [ 1 ] to express P(a\T esc ) (see Appendix [A). In Fig. 5.5 \, we confront 



in good agreement the numerically computed P(a\T esc ) and the generalized arc- 
sine law ( |A.l| l. Finally, computing numerically the integral in equation (5.5 i and 



using the relation (5.4 1 we obtain the results shown in Fig |5.5B . The numerically 



computed effective rate y is reasonably well described by the approximation (5.4 1 
up to ju/A < 10 (see also inset), while the approximation breaks down for larger 
values of fj. . 

When jU is too large, the approximation might indeed fail for two reasons. 
First, we should notice that r esc is a random quantity defined as the first time when 
the separation reach es L from a. This can be taken into account by approximating 
P{T & sc) in equation \5.5\ with an inverse Gaussian, which also comes from con- 
sidering p(f) = In (/?(/") J a) as a Brownian motion with drift. However, the gener- 
alized arcsine law used to express P(a|r esc ) looses its validity when conditioning 
on the value of the now random r esc . Second, in spite of the striking agreement in 



figure 5.5 \, we know that the approximation of p(t) as a Brownian motion with 
drift fails at short times. For example, in the case p=0, it is possible to show J2] 
12 1 that dp/d? = A coth(2p) + yAt] (t). The approximated dynamics is thus only 
reached when coth(2p) f» 1, that is at long-enough times. 




Fig. 5.5 (A) P(a\T esc ) with T esc = ln(L/a)/X(p) for different values o f pas labeled and L/a = 
2000. The solid lines correspond to the modified arcsine law prediction \A.l) obtained assuming 
a Brownian motion with drift X (p) = 2D\{\ — 2p) and diffusion coefficient A (p) /2 = P i (1 + 
2p). (B) y/y„ vs for the same values of p, the solid lines correspond to (1 — j8), see Eq. \5A\ . 
The inset is a zoom of the main figure in log-lin scale. 



22 



6 Conclusion 

In this work we have considered a dilute system of particles transported by a com- 
pressible flow. The dynamics of particles was supplemented by the binary reaction 
A + A — > with a rate jU. Particles can react only when their separation is below 
a finite interaction radius a. We focused on a situation where transport dominates 
over diffusion (high Schimdt numbers) and thus particles are only driven by the 
carrier flow. The interest of the setting of this work is that, independently of the 
initial number of particles, the standard mean-field approach always breaks down 
at sufficiently long times. Also, due to the compressible transport, strong inho- 
mogeneities appear in the particle spatial distribution, changing the local effective 
reaction rate. To describe such a system we introduced the joint «-point density 
that was shown to obey a hierarchy of transport equations. This hierarchy was then 
solved using a Lagrangian approach based on averaging over particle trajectories. 
The solutions are shown to decay exponentially at long times unlike the algebraic 
oc t~ l standard mean-field prediction. Then, it was explicitly shown that all mo- 
ments of the number of particles decrease with the same rate y. This effective rate 
depend on the statistics of the relative distance between a pair of Lagrangian tra- 
jectories. One of the main results of this work was to provide a formulation of 
the problem in terms of non-interacting Lagrangian trajectories of the flow. This 
approach allowed us to apply some standard tools of statistical mechanics and pas- 
sive (non-reacting) compressible transport to the problem of advection-reaction. 
Assuming the existence of a large deviation principle for the time that two non- 
reacting particles spend within a distance a, we expressed the effective rate in 
terms of the (Cramer) rate function. The asymptotic values of the effective rate y 
for /I — > and jj. — > °° were obtained using general properties of the rate function. 
We then introduced a phenomenological model that allowed us to understand how 
the effective rate depends on the compressibility of the carrier flow and on the 
interaction radius. Finally, we numerically studied the case of reactants advected 
by a smooth compressible Gaussian delta-correlated-in-time flow. Numerical and 
theoretical results were found to be in good agreement. 

Many aspects of this work can be easily extended to more complicated situa- 
tions. For instance we have focused only on tracers, that is particles that just follow 
a prescribed flow and whose ^-dimensional position space dynamics is given by 
equation \2.l) . This choice was made in order to simplify the presentation. Most of 
the results presented in this work remain valid when considering particles obeying 
a Newton equation, that is whose dynamics is in the 2ii-dimensional position- 
velocity phase space. For instance, the expression ( |3.8| l for the mean number of 
particle in terms of relative distances remains formally the same, but one has to 
adopt a gas kinetics approach and perform an additional integration over particle 
velocities. Also, we considered particles reacting within a sharp interaction radius. 
This hypothesis can be easily extended replacing the Heaviside function by a more 
general kernel, which could possibly depend on velocity differences or any other 
intrinsic property of the particles. 

Despite these straightforward extensions, many questions are still open. For 
instance, we showed that factorial moments of the number of particles exponen- 
tially decrease with a rate that is an increasing function of their order. However, 
only in the zero-dimensional case an explicit dependence was given (namely y„ = 
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\xn(n — \)/2j2). In d > 1, these effective rates depend on the joint statistics of sev- 
eral particle pairs. Multi-time correlations of non-independent pairs of Lagrangian 
trajectories directly intervene making difficult to evaluate the effective rates. De- 
veloping such theoretical tools is important not only for the reaction-advection 
problem but also for general questions associated with particle transport. Another 
issue that was not addressed in this work, corresponds to the case when the ve- 
locity flow is not smooth enough to guarantee the existence and uniqueness of the 
Lagrangian trajectories. This is for example, the case of turbulent flows when the 
interaction radius is greater than the Kolmogorov length scale. In such cases, even 
if the solution in terms of tracers could be in principle extended by using the sta- 
tistical description of the so-called generalized Lagrangian flow lH6l . many of the 
theoretical considerations can drastically change. For instance, in these situations 
the relative dispersion of particles is very different from the Brownian-motion- 
with-drift lognormal description used in this work. 

Another topic, not covered by this work, concerns events related to the right- 
hand side of the rate function, which corresponds to the rare situations where 
particles stay close for times much larger than the average. This is relevant for 
instance when considering the probability that particles react very fast, a situation 
that can lead to underestimating the amount of reactants needed in some practical 
applications. 

Finally, a natural extension of this work is to consider more general reactions. 
In particular, the case of coalescence/coagulation in dilute media presents many 
relevant applications to different fields, such as astrophysics or geosciences. Note 
that an alternative approach to that given by the Lagrangian trajectories consists in 
trying to close the hierarchy of transport equations by performing an appropriate 
coarse-graining. However, because of the finiteness of the interaction radius this 
approach is not straightforward. 

A Generalized arcsine law 

In this brief Appendix we summarize the generalized arcsine law derived in Ref.|lj for the 
probability density function P(a\T) of the fraction of time a spent in the negative semiaxis 
during a time interval T by a one-dimensional Brownian motion with drift 8 and diffusion 
coefficient D. Denoting Y = 8 /\/2D, such a probability reads 

P(<*\T) = \ 

with %{x) = jerfc [xj^fl^ . Notice that for zero drift T = the above expression recovers the 
standard arcsine law P(a\T) = ^/a(\ — a]. 
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